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Abstract 



We present analytical investigations of a multiplicative stochastic process that 
models a simple investor dynamics in a random environment. The dynamics of the 
investor's budget, x(t), depends on the stochasticity of the return on investment, 
£> ' r (^)> f° r which different model assumptions are discussed. The fat-tail distribution 

^ ■ of the budget is investigated and compared with theoretical predictions. We are 

■ mainly interested in the most probable value x mp of the budget that reaches a 

constant value over time. Based on an analytical investigation of the dynamics, we 
are able to predict x^pf • We find a scaling law that relates the most probable value 
fx*. \ to the characteristic parameters describing the stochastic process. Our analytical 

results are confirmed by stochastic computer simulations that show a very good 
agreement with the predictions. 

Keywords: multiplicative stochastic processes; scaling laws; investment dynamics; PACS: 
05.40.-a, 05.10.Gg 



1 Introduction 



Multiplicative stochastic processes denote type of dynamics where a variable x(t) changes 
its value in time due to a stochastic term A(t): 



x(t + 1) = \(t) x{t) (1) 
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A(t) may describe different stochastic processes, such as Gaussian white noise, uniform 
random distributions, GARCH or ARCH processes, which are explained later in this 
paper. Despite its great simplicity, the dynamics of eqn. ([1]) gained much attention in 
different fields of applications. It was, for example, used already in 1931 by Gibrat to de- 
scribe the annual growth of companies - an idea extended in different works by economists 
10, 42] and econophysicists [l|. Several theoretical aspects of stochastic processes with 



multiplicative noise have been focus of recent research in physics 20|, |3jJ, |40|, |4]J, |44j] • More- 



over, by taking into account additional couplings between these processes, the approach 
has been extended towards a generalized Lotka-Volterra model 22j, |39J which is referred 
to later in this paper. 

The importance of multiplicative stochastic processes was also reflected by several math- 
ematical investigations 12|, |l3j, |46(. In the so-called Kesten process 16J - from which 
we depart in this paper - the dynamics of eqn. (CQ) was extended by a second additive 
stochastic term a(t) independent of A(£) 



x(t + l) = X(t)x{t) +a(t) 



(2) 



This extension has the nice feature that the stochastic process is repelled from zero, 
provided some constraints on a(t) are satisfied. We will come back to this later in our 
paper. At this point, we rather want to provide some arguments why we became interested 
in this topic and why we think that this is relevant. 

Our investigations started with the question of how much "intelligence" is needed for an 
agent to survive in a noisy environment (see also 0). The example at hand is explained in 
Sect. 12.11 agents with a personal budget, x(t), participate in a simple investment scenario 
and face a return on their investment, r(t). Given some uncertainty of r(t), they have the 
choice to adjust individually the portion of the budget they wish to invest, q(t), which 
can be called their risk propensity. The strategy to choose q(t) may of course depend on 
whether the agents are able to obtain some information about the expected return on 
investment, r(t). Here, different levels of internal complexity, or "intelligence", of agents 
come into play, namely their capability to observe and to store a history of previous returns 
r(t — t), to calculate different measures from such a time series (such as trends, moving 
least squares, etc.), to detect or to forecast certain periodicities in the signal received (such 
as cycles in the market) (2s|, 27]. At this point, a wealth of different models, assumptions, 
suggestions, speculations etc. about the agent's behavior comes into play from various 
fields (financial theon 



„ 13, 12J, 130 



economics 



13, [2l|, |29j, |45[ , behavioral sciences 



43], physics 0,0,0, computer sciences 0,[3,[l8|). All these rather complex assumtions 
eventually lead to their specific outcome of the game, and can hardly be compared. So, 
before embarking on some more refined modeling assumptions about agents behavior, we 
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had to ask ourselves what could be the reference case in this simple investment scenario, 
to which later all advanced simulations can be compared. In other words: what would 
be the dynamics of that investment process without all these "intelligent" assumptions? 
It turns out that this baseline case is exactly given by the dynamics of a multiplicative 
stochastic process with an additive term as described by eqn. (T5]). Instead of a rather 
complex strategy q(t), we simply assume go = const, as the reference case, instead of an 
unknown additional influx a(t), we take a constant, but positive "income" a, and instead 
of specific economic assumptions about market fluctuations we choose r(t) from four 
different stochastic processes, which are simple, but analytically tractable. This way, in 
Sect. 12.11 we arrive at the basic dynamics of eqn. (1121) . which is equivalent to our starting 
equation 

The aim of our paper is (i) to elucidate the dynamics of eqn. (T5]) for some specific settings of 
X(t) and a(t) by means of some computer simulations, and (ii) to investigate its stationary 
properties by means of analytical treatment. This will reveal some scaling between the 
most probable value of x(t) and the parameters describing the stochastic processes. Before 
we continue in this direction, we want to shortly summarize some previous theoretical 
investigations of eqn. (|2]), in order to refer to it later. 

The dynamics of eqn. © was treated by Sornette and Cont as 

x{t + l)-x{t) _ a(t) 

xjt) "^ + A(t)_1 (3) 

If the finite difference (x(t + 1) — x(t))/x(t) is approximated by (dlogx/dt), the following 
overdamped Langevin equation for w = logx can be obtained: 

^ = a(t)e- w - \v\ + V (t) (4) 

with: 

v = <A)-l^<logA> W) = <A 2 > - (A) 2 (5) 

The first term on the r.h.s. of eqn. (j4]) describes an effective repulsion of x from zero, while 
the second term, — \v |, describes the drift towards zero. The third term, t](t), expresses 
the stochastic influences resulting from \(t). 
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Going over from the single stochastic realizations of u(t) to the probability density P(u, t), 
it was shown in 41] that the following Fokker-Planck equations can be derived 



dP(w,t) 
di 



a(t) e- w P(w,t) 
+ «log(A) 2 )-(logA) 2 ) 



((log A) + a(t) e~ w ) 
ld 2 P( Wl t) 



dP{w,t) 
dw 



(6) 



In accordance with eqn. 



2 dw 2 

the first term of eqn. (jSJ) describes the decay on w, the second 
term indicates the drift of the process and the third is a diffusion term with the diffusion 
constant: 

D = <(logA) 2 )-(logA> 2 (7) 

We recall that without the additive stochastic force, i.e. a(t) = 0, eqn. dSj) results in 
a simple Fokker-Planck equation for x(i), related to the stochastic eqn. (pQ), with the 



log-normal distribution as limit distribution 3l|, |41 

11 f 1 



P(x(t)) 



VnDt x(t) 



cxp 



Dt 



[logx(t)-(logA> t] 



Considering, however, an additive stochastic force a(t) ^ 0, it was already noted in 16 



that for large x instead of a log-normal distribution now a power-law distribution results, 
provided that (log A) < 0: 

P(x) oc x- (1+At) (9) 
The exponent fi satisfies the conditions 

(A"> = 1; fi 



2 I (log A) | 



((logA)2)-(logAr 



2 |u| 



(10) 



It was shown 4lj that such power-laws result under quite general assumptions about 
multiplicative stochastic processes with repulsion at the origin which generalizesprevious 



;s_p 

results by [20j. We note that an economically motivated entry/exit dynamics [3, |32|, [36 
as a boundary condition of the multiplicative stochastic process also leads to power dis- 
tributions. 



In the more general case of eqn. (jHJ), we will in Sect. 13. II provide a stationary solution for 
P(oj) at least for the case of a constant a. We will further use this solution to derive a 
result for the most probable value of the multiplicative stochastic process. In order to test 
the analytical predictions, we will compare them with stochastic computer simulations, 
which are described in the following section. 



1 Note that this equation is different from the original one in 4l| by a factor 2 in the diffusion term 
which later affects the definition of fi in eqn. (jfOp 
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2 Computer simulations for different random envi- 
ronments 

2.1 Investment model 

Instead of solving the stochastic differential equation ([6]) numerically, in this section we 
focus on the individual realizations of the stochastic process as described by eqn. (EJ. 
Hence, we have to specify the stochastic processes behind the variables \(t) and a(t), 
which is done in the following. As it became clear from the previous discussion, the 
additive stochastic term acts as a repulsion of the dynamics from zero, i.e. without a(t) 
after some time every stochastic realization of x(t) will reach the value zero, i.e., the 
exit for the process. So, the meaning of a(t) is simply to keep the dynamics "alive" by 
preventing it from reaching zero. To simplify the dynamics, instead of a time-dependent 
value a(t) we simply choose a small, but constant positive value of a, which results in the 
same effect. 

For the multiplicative stochastic term the following assumption was made: 

X(t) = l+r(t)q(t) (11) 

In order to give the stochastic dynamics of eqs. ([2]), ffTT]) some interpretation, let us 
consider the following: x(t) shall represent the budget (wealth, liquidit y) o f an agent, who 



intents to invest some portion < q(t) < 1 of its budget into a market |26|,|27|. Depending 
on the market performance, this investment may result in a loss or a gain. Hence, there 
exists a return on investment, Rol, r(t), which is either positive or negative. It should be 
noted that the Rol has to satisfy r(t) > — 1 as a lower boundary, as investors can never 
loose more than what they have invested. Because there is no limit for potential profits, 
in principle there is no upper boundary for r(t). 

In a real investment scenario, the Rol can be determined from a time series of a real 
market that is influenced by several thousands bids and asks. This would involve to 
model the market dynamics explicitely. In the spirit of other investigations in economics 
and econophysics [l| we have chosen instead to model the market return r{t) by means 
of different stochastic processes, i.e. to keep its dynamics independent of the investment 
of the agent. Further, we have bound the Rol to values between -1 and +1. The upper 
boundary was chosen to obtain a mean of zero for r(t), which allows to better understand 
the basic dynamics. The following distributions for r(t) are discussed: 

1. B{ — 1, 1}, a binary stochastic switch between only two states, -1 and +1 
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2. U(—l, 1), a uniform distribution, where every possible value between -1 and +1 has 
the same probability to be chosen 

3. N(0, 0.1), a normal distribution with mean and standard deviation 0.1, i.e. values 
close to the mean are more likely to be picked 

4. ARCH(l), a distribution, where possible values between -1 and +1 are chosen by 
considering some correlations between times t and t — 1 (autoregressive conditional 
heteroskedasticity) 0,0]. 

These stochastic processes were not chosen in the first place because they capture real 
economic processes e.g. on financial markets, but because they later allow for analytical 
calculations of some properties of the probability distributions, as shown in Sect. 13.21 
However, we note that the four types cover different degrees of complexity in the stochastic 
process, as their spectrum of values (discrete, continuous) and their range of values (rather 
broad for the uniform distribution, but narrowed and centered for the normal distribution) 
vary and even correlations between different time steps are taken into account (ARCH(l)). 

The investment dynamics is then described by the multiplicative stochastic process 

x(t + 1) = x{t) [1 + r{t) q{t)\ + a (12) 

where the positive constant a prevents the investor from going bankrupt. 

The challenge for an agent in this very simple investment model is then to adjust its 
portion of the budget to be invested, q(t), given some (observed) market returns r(t). 
Low values of q refer to a risk averse investment strategies, while higher values may be 



more risky, but also more rewarding with respect to the budget. In [27J we have discussed 
different scenarios for agents to adjust their risk propensity, q(t). In this paper we focus 
more on the stochastic dynamics and therefore just choose q(t) = q^, which is a small but 
constant value. This can serve as a reference case for more complex investment strategies 



27 



At this point, we wish to note that the dynamics proposed in eqn. (fT2l is much simpler 



than the comparable starting equation in the generalized Lotka-Volterra model 22|, 139 
mentioned in Sect. (TJ 

Xiit + 1) = Xi(t) [1 + \{t)} + ax{t) - Cxiit) (13) 

which assumes additional couplings between different individual stochastic processes, 
Xi(t). Instead of a small, but constant income a the term ax(t) considers a global cou- 
pling via the mean budget x(t) of all agents (which may be related to general publicly 
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funded services). Moreover, the third term Cxi(t) describes direct interactions between 
different agents, as C is a function dependent on other Xj. These may account for com- 
petition for limited resources and saturation effects in the dynamics. Even if, after some 



approximations discussed in |22[, these additional influences may be small, they still af- 
fect the general solutions for the underlying probability distributions as we will discuss in 
Sect. 13 . 1 L If we recast the differences between eqs. ( [121) . ( |T3l) using the general framework 



of multiplicative processes [33 



Ax{t) = 7]{t)G[x{t)} + F[x(t)\ (14) 

where rj(t) is a stochastic variable, then we have for our model, eqn. (1121) G(x) = x and 
F(x) = a, whereas for eqn. (fl3|) G(x) = x and F(x) = a(l — x) results |22j. This will 



eventually affect the specific exponents of the stationary probability distributions, P s (x). 
2.2 Simulation results 

Here, we present stochastic computer simulations of eqn. (Tl2|) for different distributions of 
r(t) as described in the previous section. Initially, x(0) = 10 holds for the agent's budget, 
go and a are kept constant during each simulation. The distributions were realized using 
lO'OOO agents. In the case of the time evolution of the most probable value x mp , the plotted 
values were calculated averaging over the x mp obtained for distributions resulting from 10 
different realizations of the simulation. 

Figure [JJ shows the distribution of the investor's budget at different time steps. Starting 
from the delta distribution at x(0) = 10, the distribution disperses over time until it 
reaches a stationary state. The number of time steps needed to reach this stationary state 
depends on the initial conditions, the distribution of r(t) and the additive constant a. For 
the conditions in Figure [JJ this takes about 10 4 iterations. 

One can clearly see that the stationary distribution is characterized by a fat tail described 
by a power-law distribution, as shown in Figures [21 This is in accordance with previous 



investigations 20|,|40|] and was already discussed in Sect. [TJ We have determined the scaling 
exponent \x of the power law, eqn. [HI from the simulation data for different stochastic 
processes and will later compare it with our analytical investigations. 

In the following, we are more interested in the clear maximum of the stationary distribu- 
tion which gives the most probable value, x mp . Figured (right) shows the evolution of x mp 
over time for three different additive terms. The most probable value increases with time 
until a balanced state is reached, where the dissipation described by the negative drift 
towards zero and the constant influx a compensate on average. It can be noticed that the 
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Figure 1: Investor's budget probability distribution P(x, t) for different time steps, assum- 
ing x(0) = 10 and r(t) € B{ — 1, 1}, q = 0.1 and a = 0.5. The probability was estimated 
from N = 10 4 runs. 




Figure 2: Investor's budget stationary probability distribution P s (x) (estimated from fre- 
quencies after t = 10 4 ): (left) Binary stochastic return distribution r(t) = B{— 1, 1}. 
(right) Uniform stochastic return distributions r(t) = U( — l, 1). In both cases, x(0) = 10, 
q(t) =0.1 and a = 1. Data is binned in logarithmic intervals of the same size. The dashed 
lines show the theoretical prediction of these curves by Eqn. (l2"Tj) . derived later. 



standard deviation of x mp increases with a. Larger a lead to larger values for the budgets, 
this in turn leads to an amplification of the fluctuations due to the multiplicative nature 
of the process. Figure [3] shows the evolution of x mp for three different values of the risk 
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propensity g . Obviously, the larger q, the smaller x 
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Figure 3: Most probable budget value vs. time: (left) for q = 0.1 and 3 different addi- 
tive constants a, (right) for a = 0.5 and 3 different constant risk propensities go- Other 
parameters and settings as in Figure E^right). 



Figure [3] suggest that there is a scaling between the most probable value x mp and the 
parameters characterizing the stochastic dynamics, eqn. (1121) . namely go and a. This scal- 
ing is investigated numerically in the following and will later be confirmed by analytical 
investigations. Figures HI [5] show (x mp ) at a fixed time, t = 10 4 , for the four different 
realizations of the stochastic process, r(t), discussed in Sect. 12.11 In all four cases, the 
results, plotted against the variable ajq\ clearly show a straight line, which allows for the 
scaling: 



x 



mp 



(15) 



By means of analytical investigations of the stochastic multiplicative dynamics, we now 
want to confirm the findings of the numerical experiments. 



3 Analytical investigations 
3.1 Stationary state 

In order to derive the observed scaling law, eqn. (j!5p . we start with the Fokker-Planck 
eqn. (JB])o Using the approximation of a constant value of a and the notation of eqn. (JTj) 



J In the appendix, we point to the direct treatment of the stochastic dynamics, eqn. (|2|). 
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Figure 4: Most probable budget value vs. scaled variable a/ql for the binary stochastic 
process r(t) G {—1,1} (o) and for the uniform stochastic process r(t) G U(—l,l) (A). 
Each set was plotted by varying a and q® over the range of [0.1,0.9] in 0.1 increments, 
giving a total of 81 data points per combination, , averaged over lOO'OOO Monte-Carlo 
simulations and 10 runs. The value of the slope found for the numerical simulations is 
for the binary stochastic process c = 0.961 ± 0.006 and the uniform stochastic process 
c = 3.149 ±0.014 




Figure 5: Most probable budget value vs. scaled variable a/<7o : (left) for the normal 
distributed stochastic process r(t) G iV(0, 0.1). The value found for the slope of the 
numerical simulations in this process was: c = 102.17 ± 0.701, (right) for the ARCH(l) 
process (a — a i — 0.1), with a slope of c = 9.089 ± 0.029 
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for the diffusion constant D, we find for the stationary solution of eqn. 



-<9„ 



({\og\) + ae- w )P s 



w 



^-d w P s (w) 



This results in the stationary solution: 



P s (w) = TV exp 



2 (log A) w — 2ae 
D 



(16) 



(17) 



with normalisation A/". Using w = log a;, the corresponding stationary probability distri- 
bution P s (x) is recovered by the chain rule: 



P 3 (x) = P s (w(x)) — = Afx r> 1 

ax 



The normalization is explicitely calculated via: 

2a 



exp 



2a 



Dx 



2(logA; 



TV / 
'o 



exp 



dx = M[^- 
D x I V 2a 



2(log A) 
D 



2 (log A) 
D 



(19) 



where T(y) is the Gamma function. Redefining 

2 (log A) 



D 



(20) 



we eventually find the stationary distribution in the normalized form: 



r{ii) 



PJx) 



x exp 



2a 



Dx 



(21) 



A graphical visualization of this function can be seen in Figs. |5] together with the results 
obtained from the averaging of the stochastic computer simulations. It is worth noticing 
that the Gamma function diverges for /x — > 0. Hence, the condition of (log A) < is 
indeed important for the existence of the stationary solution. One should also notice that 
for large x eqn. (12TI) reduces to the power-law distribution, eqn. (jnj). 

We note that eqn. (121]) is a special form of the general solution 



Ps(x) 



G*(x) eXP V 



F{x') 
G^x 7 



dx' 



'221 



obtained in 22|, [33(|. If we use F(x) = a and G(x) = x in accordance with the stochastic 
process defined in eqn. ( 1T2|) . this would lead to the (non-normalized) solution P s (x) = 
x~ 2 exp(— 2a/ Dx). This is in agreement with eqn. (1211 because of fi ~ 1 in our case. The 
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solution discussed in 22J is however different from ours because F(x) = a(l — x) and 
G(x) = x was used there, which eventually lead to the stationary distribution P s (x) = 

x -2(l-a/D) exp ( 

—2a/Dx), and consequently tou = l + (2a/D). Only in the limit of a <C D 
(which may hold for the case discussed in 22| because of e.g. a ~ D/4, but is hardly 
satisfied in our model, because for example a = 0.5, D = 0.01 in Fig. CD), both solutions 
tend to converge. Apart from these differences, the emergence of the stable power-law 



distribution, eqn. (jHJ), for large x was already discussed in [3, 20, 33, il, 37, 38, 39, 40 



11 



3.2 Scaling of the most probable value 

In the following, we are less interested in the power-law behavior of eqn. (1211) . but mainly 
in the most probable value of the process, x mp , which correspond to the peaks of the 
distribution shown in Figs. [21 From the extremum condition 



= d x p s (x) 

we find with eqn. (|21|) for the most probable value: 

a 



x 



mp 



D-(logA) 



(23) 



(24) 



(Note again the difference to the result x mp = a/(D + a) obtained in [22j - for a different 
stochastic process). 

In order to derive the scaling of x mp on the parameters a and go, we have to calculate the 
diffusion constant D, eqn. (JTj), and, hence, have to specify the stochastic process, X(t), 
eqn. (ITT1) . or r(t) respectively. From eqs. (J7J), ( fill we find with q(t) = qo 

v2 _ 2 II 2\ / \2N 



D = <log 2 (l + g r)) - (log(l + q r)f w q 2 «r 2 ) - (r) 2 ) 



(25) 



Because of (r) = 0, this results in first approximation in the following expression for the 
most probable value: 



Note that this yields a good approximation only for (logX) ~ and small values of go- 

For the specification of r(t) we refer to the four different distributions listed in Sect. 12.11 
and already used in the computer simulations of Sect. 12.21 We find the following results: 

1. If the returns, r(t), are randomly drawn from the binomial distribution B{ — 1; 1}, 
(r 2 ) = 1 holds and the scaling is obtained as: 



x 



mP <ll 



(27) 
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2. If the returns are randomly drawn from the uniform distribution U(— 1, 1), we find 

(r 2 ) = \ J J 2 dr = i (28) 

and for the scaling 

x m p 3—^- (29) 
% 

3. If the returns are randomly drawn from the Gaussian distribution iV(0,0.1), the 
deviation is given by (r 2 ) = cr 2 = 0.01 and the scaling is obtained as 

x mp = 100^ (30) 
% 

4. If the returns are randomly drawn from the ARCH(l) process, we find with a = 



ai = 0.1 [23j, P- 78] 

(r 2 ) = a 2 = = 0.111 (31) 
1 — OLi 

which leads to a scaling of 

x m p 9— „ (32) 

We note that the nature of the stochastic process in all four cases is reflected only in 
the different prefactors c, while the scaling function x mp = ca/q^ remains the same. 
Comparing the analytical results with the computer simulations presented in Sect. 12.21 
one can see that both the scaling function and the values of the prefactors c are in perfect 
agreement. This shows indeed that the approximations made for the analytical treatment 
were appropriate. 



4 Conclusions 

In this paper, we have derived an analytical expression for the stationary probability dis- 
tribution of an investor's budget x, eqn. (1211) . when investing in a random environment. 
The nature of the underlying stochastic process was considered in four different distri- 
butions for the return on investment (Rol), r(t). Assuming that in every time step the 
investor invest a constant portion q of his current budget on which he receives an Rol 
and further receives a very small amount a as a constant income, we have shown that the 
most probable value of the investor's budget scales with the forementioned variables as 



x 



mp 



a/(qo (r 2 )), eqn. ( 1261) . This result was confirmed both by analytical investigations 



and extensive computer simulations for the four different stochastic processes. 
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At the end, we want to comment on the range of validity for the scaling obtained. One of 
the underlying assumptions of our model was a constraint of the Rol to values between 
(— 1, +1), where -1 means a complete loss of the investment and thus is a reasonable lower 
bound, whereas +1 means a doubling of the investment, chosen for reasons of better 
tractability (i.e. (r) = holds, for example). These constraints for r(t) used during the 
computer simulations may indeed result in deviations from the theoretical prediction, 
eqn. (EE}, if the underlying stochastic process for r(t) frequently gives values outside the 
interval ( — 1,+1), which then need to be discarded. This can be the case both for the 
normal and for the ARCH distribution. 

For our computer simulations, we have thus chosen small values of a 2 and ao, cm, respec- 
tively, to control the width of the distribution and the "outliers". In fact, the larger these 
values, the less is the agreement between the "truncated" computer simulations and the 
theoretical approximation based on the full range of r values. 

Nevertheless, this argument does not restrict the value of the scaling obtained in eqn. 
(1261) . which is still valid. In order to deal with broader distributions for the Rol, one has 
two options: 

(i) The computer simulations can be repeated on a larger interval for r(t) G (— C, +C) 
which controls the number of "outliers" . This will however not change the principal 
insights derived in this paper. 

(ii) The analytical prediction can be improved by dealing with truncated distribu- 
tions. This basically affects the calculation of (r 2 ). Assuming the constraint r(t) G 



-1, +1), one finds for example for the truncated normal distribution 35 



which holds also for larger a 2 . For truncated ARCH or GARCH distributions, the 
situations are more complicated, here we refer to the literature 



Eventually, we wish to comment on the investor dynamics proposed in this paper. These 
were by purpose related to multiplicative stochastic processes, to allow for analytical 
insights and conclusions about the influence of different stochastic distributions. Despite 
its simplicity, the dynamics of eqn. ( 1121) is not restricted to artificial scenarios. In fact, the 
dynamics of the Rol, r(t), can be taken from real time series, instead of being modeled 
by a stochastic process. The most challenging application, however, is in the dynamics 
of the variable q(t) which decides about the portion of the budget to be invested. While 
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risk averse agents may tend to lower q, risk seeking strategies may go for higher q and 
for higher yields in the lucky case. Throughout this paper, q was set to a constant but 
small value, independent of individual characteristics. Any realistic investment scenario 
deals with the question to adjust the risk propensity q(t) in time based on the observation 
of previous r(t). For the model given, an artificial intelligence approach to this question 
is presented in 22J. As a next step, the investment dynamics can be extended towards 
a portfolio scenario, where both r and q become multidimensional variables, to allow 
different investment strategies for different assets. 
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Appendix 



Instead of dealing with the probability distribution P(x) as in Sect. [31 one can try to 
treat the stochastic dynamics of eqn. ([2]) d irectly. For this case, we can present at least a 
formal solution using the ^-transform |l4j. Rewriting eqn. (T5]) as 

a = x(t + 1) - Xx(t) 



the ^-transform leads to 

oo ^ 
n=0 

where X(z) is given by 



z - 1 



z [X(z) - z(0)] - AX(z) 



X( ^ ) = x(0) ^A +a (,-!)(, -A) 



Using the inverse tranform 

X(^" 1 ) = z(0) 
the solution for can be found as 

X(t) = 



1 - A 



a 1 



1 - Xz 1 - A 1-2 



A f 



x(0) 



A*x(0) + a 



1-A 
1-A* 

1-A 

t-i 



1-A 



A'x(0) + a J2 X * 



(34) 



(35) 



(36) 



(37) 



(38) 
(39) 
(40) 
(41) 



From this solution, we see that the decisive condition on A for a well-defined solution 
reads: 

|A| < 1 (42) 

which agrees with the finding of (log |A|) < obtained from the treatment of the proba- 
bility distribution P(x). 
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